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We consider the conductance of an Andreev interferome- 
ter, i.e., a hybrid structure where a dissipative current flows 
through a mesoscopic normal (N) sample in contact with two 
superconducting (S) "mirrors". Giant conductance oscilla- 
tions are predicted if the superconducting phase difference <j> 
is varied. Conductance maxima appear when <j> is on odd 
multiple of 7r due to a bunching at the Fermi energy of quasi- 
particle energy levels formed by Andreev reflections at the 
N-S boundaries. For a ballistic normal sample the oscilla- 
tion amplitude is giant and proportional to the number of 
open transverse modes. We estimate using both analytical 
and numerical methods how scattering and mode mixing — 
which tend to lift the level degeneracy at the Fermi energy 
— effect the giant oscillations. These are shown to survive 
in a diffusive sample at temperatures much smaller than the 
Thouless temperature provided there are potential barriers 
between the sample and the normal electron reservoirs. Our 
results are in good agreement with previous work on conduc- 
tance oscillations of diffusive samples, which we propose can 
be understood in terms of a Feynman path integral descrip- 
tion of quasiparticle trajectories. 



I. INTRODUCTION 

Recently considerable attention has been devoted to 
mesoscopic superconductivity, i.e. to the transport prop- 
erties of mesoscopic systems with mixed normal (N) and 
superconducting (S) elements, where new quantum inter- 
ference effects have been discovered. Novel physics ap- 
pear in such systems because electrons undergo Andreev 
refiectionsEl at the N-S boundaries, whereby the macro- 
scopic phase of a superconducting condensate is imposed 
on the quasi-particle wavefunctions in the normal regions. 
If transport in the normally conducting part of the sam- 
ple is phase coherent, there is a possibility that interfer- 
ence between Andreev scattering at two (or more) N-S 
interfaces makes the conductance of the hybrid system 
sensitive to the phase difference cf> between the super- 
conducting elements; in this case one may describe the 
system as an Andreev interferometer. 

This paper is concerned with a theoretical description 
of hybrid mesoscopic systems of the Andreev interferom- 
eter type. In particular we are interested in the nor- 
mal conductance as a function of the phase difference 
between the condensates of two separate superconduct- 
ing elements acting as "mirrors" by reflecting the qausi- 
particles in the normally conducting element, which in its 



turn is connected to two electron reservoirs as schemat- 
ically shown in Fig. |l|. The normal electron transport 
may be in the ballistic regime or in the diffusive regime; 
both cases will be discussed. In addition we will make the 
important distinction between the cases when potential 
barriers or (sharp) geometrical features serve as "beam 
splitters" at the junctions between the leads and the nor- 
mal element and when the passage between leads and 
sample is unhindered by quantum-mechanical scattering 
between distinct quasi-classical trajectories at these junc- 
tions. 




FIG. 1. Schematic picture of an Andreev interferometer 
consisting of a normal (N) metal (diffusive transport regime) 
in contact with two superconducting elements (S), which are 
characterized by the phases cj>i and 02 of their respective order 
parameter. The normal metal are in contact with two reser- 
voirs Ri and R2 via leads. The thick lines at the junctions 
between the leads and the normal metal represent potential 
barriers, which act as beam splitters partially reflecting quasi- 
particles impinging on the junctions, cf. Fig. 2. If transport 
is phase coherent quasiparticles at the Fermi level (zero exci- 
tation energy) are phase conjugated after Andreev reflection 
at the N-S interfaces, so that Andreev reflected holes (dashed 
line) retrace the path of the incoming electrons (full line) and 
vice versa. 

The rest of this introduction will be divided into two 
parts: (i) a general introduction to the subject of An- 
dreev interferometry and (ii) a qualitative discussion 
based on quasi-particle trajectories which makes it pos- 
sible to understand the main features of the conductance 
oscillations of various types of Andreev interferometers a 
a function of the superconducting phase difference. 
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A. Origin of conductance oscillations in Andreev 
interferometers 

Already in the early 80's Spivak and Khmel'nitskii 
showeda the weak localization corrections to the con- 
ductance of a diffusive sample containing two supercon- 
ducting mirrors to be sensitive to the superconducting 
phase difference. The effect can be understood in terms 
of the usual interpretation of weak localization as due 
to coherent backscattering. The interference of probabil- 
ity amplitudes for classical quasi-particle trajectories (or 
"Feynman paths" ) bouncing off both mirrors will depend 
on the phases of the respective condensates. Considering 
a closed diffusive path touching both S-N interfaces - 
where electrons will be reflected as holes and vice versa 
— the interference between quasi-particles moving in op- 
posite directions, clockwise and counterclockwise around 
the path, results in a phase difference of 2<j) between the 
interfering amplitudes, i.e. twice the phase difference 
between the two superconductors. This is because the 
phase picked up due to Andreev reflections off the two 
mirrors is ±0 depending on whether the motion is clock- 
wise or anticlockwise. It follows that the weak localiza- 
tion correction to the conductance of a normal sample 
with two superconducting mirrors has a component that 
oscillates with a period equal to tt as the phase difference 
between the superconductors is varied. 

In the beginning of the 90's, a dependence on the 
phase difference <p was discovered not only for the conduc- 
tance fluctuations but for the main conductance as well. 
Not only conductance fluctuations but the ensemble- 
averaged conductance itself can therefore be controlled 
by the rPUpse difference between two superconducting 
mirrors.B"l3 Hybrid S-N systems (Andreev interferome- 
ters) which show such a behaviour at very low temper- 
atures have lead-sample junctions which act as "beam 
splitters" in the sense that a quasi-particle approaching 
the junction along a quasi-classical trajectory is only par- 
tially transmitted. Hence, a beam splitting junction has 
the effect of partly reflecting a quasi-particle coming from 
one N-S interface towards the second N-S interface, as il- 
lustrated in Fig. This has the important consequence 
that when a quasi-particle finally leaves the sample to 
contribute to the current there is a certain probability for 
it to have interacted with both superconducting mirrors. 
To be specific, an electron entering the sample from one 
reservoir, R\ say (referring to Fig. [l]), may follow a tra- 
jectory (full line) where first it is reflected as a hole by one 
mirror, say S ( f >1 , then it returns (dashed line in Fig. [y) to 
bounce off the same junction through which it entered, 
gets reflected towards the second superconducting mir- 
ror Sfo where it is Andreev reflected as an electron, and 
finally it passes (full line) through the junction to the 
second reservoir i?2 now carrying information about the 
difference <p — 4>\ — (f>2 between the phases of the two 
superconducting mirrors (the difference appears because 
the phase picked up on Andreev reflection differs in sign 



between an electron and a hole). It followsoB that such 
trajectories contribute a term to the conductance that 
oscillates with period 2tt (rather than n) as function of 
<t>. 




FIG. 2. A potential barrier at the junction between the 
normal metal and the lead to a reservoir splits the quasipar- 
ticle beam coming from one of the superconducting mirrors 
(cf. arrows). This makes it possible for quasiparticles having 
undergone Andreev reflection at both mirrors to contribute to 
the current even if their excitation energy is zero and therefore 
reflected hole (electron) excitations retrace the paths of the 
electron (hole) excitations. It follows that the conductance 
may depend on tpi — <f>2 even at zero temperature (see text). 

As we have indicated above the influence of the su- 
perconducting phase difference on the conductance of an 
Andreev interferometer structure is an interference ef- 
fect. The macroscopic phases of the superconducting 
condensates, or — using a different language — of the 
order parameter or of the gap function of the respective 
superconducting mirror are imposed on the microscopic 
wavefunctions of the electron- and holelike quasi-particles 
when they undergo Andreev scattering at the N-S bound- 
aries. The dominating role of these scattering phases is 
due to the effect of compensation of the phases gained 
along the electron- and hole sections of the trajectories 
connecting the Andreev reflection. Returning to Fig. [I], 
we note that it illustrates how an electron (hole) with 
energy infinitely close to but above (below) the Fermi 
energy follows a trajectory [full (dashed) line] towards 
an N-S interface. When it gets Andreev reflected as a 
hole, conservation of energy and momentum makes its 
energy infinitely close to but below (above) the Fermi 
energy and the hole (electron) retraces the path [dashed 
(full) line] of the incoming electron (hole) . In this way the 
phase acquired by an electron is "eaten up" as the hole 
retraces the electron path in the opposite direction and 
the net change of phase is due to the Andreev reflection 
only. 

The possibility of phase compensation only exists for 
quasi-particles whose energies are very close to the Fermi 
energy. Because energy and momentum is conserved in 
the Andreev scattering process a quasi-particle with ex- 
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citation energy E measured from the Fermi energy 
is reflected as a hole of energy — E and in a direction 
that differs from the incoming path by an angle of order 
Ejtp. This implies that for finite quasi-particle excita- 
tion energies the phase compensation will not be com- 
plete. Since the dominating role of the superconducting 
phase difference is lost when the uncompensated phase 
along the quasi-particle trajectory connecting the two su- 
perconducting mirrors is of order 27r it follows immedi- 
ately that only quasi-particles whose excitation energies 
are less than a critical energy E c may contribute to the <f>- 
dependent part of the conductance. For ballistic samples 
E c ~ Tivp/L, while E c ~ hD/L 2 in the diffusive trans- 
port regime, where it is known as the Thouless energy 
(vp is the Fermi velocity, D is the diffusion constant). 

The restriction on quasi-particle excitation energies 
translates into a temperature dependence, where the 
Thouless energy sets the characteristic temperature scale. 
Nazarov et ala and Volkov et al.U have for example sug- 
gested a "thermal mechanism" that gives a large ampli- 
tude of the 27r-periodic conductance oscillations with <p at 
temperatures close to the Thouless temperature. Their 
result is due to a dependence of the effective diffusion 
coefficient on the energy of the quasiparticles in a hybrid 
S-N-S sample and will be further discussed below. 

In addition to dephasing effects due to finite excita- 
tion energies phase coherence may be broken by inelastic 
scattering. The interference effects described can there- 
fore only be observed if the length L of the normally 
conducting part of the sample is at most of the order of 
the phase breaking length L$ or the normal metal co- 
herence length Lt, whichever of them is smaller. In the 
ballistic transport regime Lt — fivp jkoT , while in the 

1/2 

diffusive transport regime one has Lt = (hD/ksT) 

A large number of experimental and theoretical 
investigationaiEj followed the early work on the tunable 
conductance of mesoscopic samples of the Andreev inter- 
ferometer type. For diffusive samples the amplitude of 
the conductance oscillations has been found to be large 
in the sense that it is comparable to the conductance in 
the absence of superconducting elements. The conduc- 
tance maxima usually appear at_even multiples of tt. As 
discussed by Kadigrobov et alx3 the situation is quite 
different for ballistic Andreev interferometers, where the 
conductance oscillations may be giant — i.e. the oscilla- 
tion amplitude may be much larger than the conductance 
in absence of superconducting mirrors. The system dis- 
cussed in Rcf. [l8] is shown in Fig. |l|; the normal part of 
a hybrid S-N-S system is weakly coupled to two normal 
electron reservoirs and hence the dissipative current flows 
from one reservoir to the other via the normal metal ele- 
ment. Two low-transparency barriers form the junctions 
between sample and leads (going to the reservoirs) and 
act as beam splitters in the sense outlined above. 

The giant conductance oscillations arises because the 
structure considered in Ref. |l^ permits resonant trans- 
mission of electrons and holes via the normal part of the 



sample. Resonant transmission occurs when the spatial 
quantization of the electron-hole motion in the meso- 
scopic normal element leads to allowed energy levels co- 
inciding with the Fermi energy (at zero temperature and 
small bias voltage the energy of the electrons incident 
from the source reservoir is equal to the Fermi energy). 
It follows from the semiclassical Bohr-Sommerfeld quan- 
tization rule [cf. Eq. ( |l2|) below] that all the N± con- 
ducting transverse modes in the normally conducting el- 
ement have one quantized level at the Fermi energy if 
the phase difference <f> between the two superconductors 
is equal to an odd multiple of tt. This means that for 
<f> = Tr(2k + 1), k = 0,±1,±2... each transverse mode 
can resonantly transmit electrons, and hence N± trans- 
verse modes contribute to the resonance simultaneously. 
As a result the amplitude of the conductance oscillations 
reaches the maximal value G max = N±2e 2 /h when <f> is 
an odd multiple of n (giant oscillations). 



B. Understanding conductance oscillations in 
Andreev interferometers in terms of Feynman paths 

In all experimental and theoretical studies of Andreev 
interferometers three types of quasi-particle scattering 
mechanisms (in various combinations) have to be taken 
into account. Scattering of charge carriers can be due to: 

1. potential barriers or geometrical feautures (beam 
splitters) at the junctions between the mesoscopic 
sample and the leads to the electron reservoirs 

2. impurity scattering inside the mesoscopic region 

3. non-Andreev (normal) reflection from potential 
barriers at the N-S interfaces. 

Here we shall emphasize the crucial role played by 
beam splitters in distinguishing between different types 
of oscillation phenomena. Therefore we choose to sepa- 
rately discuss two different types of hybrid S-N-S struc- 
tures: those with- and those without beam splitters. In 
particular we will show below that the presence of beam 
splitters is necessary for conductance oscillations with <f> 
to appear in the limit of vanishing temperature. 



1. Andreev interferometers without beam splitters 

In the absence of beam splitters quasiparticles are bal- 
listically injected into the mesoscopic sample along quasi- 
classical trajectories without suffering any quantum- 
mechanical scattering between trajectories at the junc- 
tions between sample and (leads going to the) reservoirs. 
In this case the quasi-particles therefore freely pass the 
contact region without undergoing reflection. It is not 
difficult to convince oneself that in such a system there 
are no low energy quasi-particle trajectories connecting 
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the reservoir (or reservoirs) and both superconductors. 
This is because a quasi-particle with vanishingly small 
excitation energy is perfectly backscattered at the N-S 
boundaries in the sense that the angle of Andreev scat- 
tering is equal to n. Therefore such a trajectory can not 
connect more than two bodies (say, the reservoir and one 
of the N-S interfaces). Of course, if the energy increases, 
then due to inelastic Andreev reflection (an electron with 
energy E is transformed into a hole with energy — E; the 
total energy is conserved since a Cooper pair of zero en- 
ergy is created) the back-scattering is not perfect and, in 
contrast to the case when E = 0, the angle of reflection 
differs from ir by a value a ~ E/ep. An interference 
effect involving the condensate phases of both mirrors 
is now possible since an electron Andreev-reflected as a 
hole at the S-N interface follows a different trajectory 
than the impinging electron and hence has a finite prob- 
ability not to reach the injector region. In this case, as 
shown in Fig. ||, it is possible that the trajectory will 
reach the second superconductor before finally escaping 
to a reservoir. One may readily evaluate the role of the 
described inelastic Andreev reflection in the formation of 
phase sensitive trajectories for both the ballistic and the 
diffusive case. 




FIG. 3. At finite excitation energies E the path of an An- 
dreev reflected hole (electron) in a ballistic system deviates 
by a finite angle a(E) from the path of the incoming elec- 
tron (hole). If E is sufficiently large — but not otherwise — 
one quasiparticle trajectory may therefore, as shown here, hit 
both superconductors. 

In the ballistic case an injected electron which is An- 
dreev reflected as a hole will not directly return to the 
injector if the distance Lq to the superconducting mir- 
ror is large compared to the injector opening do. The 
precise criterion is that it will not return if aLg > do, 
where a E/ep (see Fig. ||). If one takes into ac- 
count that the excitation energy as explained above is 
limited by the (ballistic) Thouless energy E c in order for 
phase coherence to be maintained, one concludes that 
an interference effect involving both Andreev mirrors 



is possible only if the injector opening is smaller than 
the electron deBroglie wavelength; this is because since 
a(E) < a(E c = hvp/L ) it follows that aL < d if 
do < Xb- (For a degenerate electron gas the deBroglie 
wavelength is equal to the Fermi wavelength.) The in- 
evitable conclusion is that an interference mechanism in- 
volving thermally excited quasi-particles cannot play a 
role in realistic experiments using ballistic samples. Un- 
der these circumstances the effect of scattering by im- 
purities inside the mesoscopic sample is decisive for the 
desired interference phenomenon involving two supercon- 
ducting mirrors to occur. In other words — in the ab- 
sence of beam splitters — we need to consider a meso- 
scopic sample in the diffusive transport regime. 

In the diffusive case interference between Andreev scat- 
tering at two spatially separated superconducting mir- 
rors may occur if the mirror-reflected trajectory diverges 
from the incident trajectory by more than a de Broglie 
wavelength, Xb, which we take to be the width of any 
particular trajectory. In this case we may say that by 
inelastic Andreev reflection the reflected quasiparticle is 
sent into a different, classically distinguishable trajectory. 
When the separation becomes greater than Xb this tra- 
jectory interacts with a different set of impurities which 
will take the reflected particle on a diffusive random walk 
along a completely different Feynman path. As the dis- 
tribution of trajectories is homogeneous in the diffusive 
limit, there is a finite probability for the trajectory (which 
starts from a reservoir) to include points with Andreev 
reflections from both superconductors. This implies that 
the criterion for the incident and reflected trajectories to 
be sufficiently separated after a diffusing length of Lo 
is ahp, > Xb- Since the angle a ~ E/ep this can be 
converted to a criterion for the excitation energy of the 
form E > E c . We recall that interference is destroyed 
for energies E ^> E c . Hence we conclude that there is 
a distinct group of quasiparticles with excitation energy 
around the Thouless energy ~ E c for which there is an in- 
terference effect controlled by the superconducting phase 
difference </>. As a result the temperature dependence 
of the conductance oscillations is non-monotonic. The 
amplitude of the oscillations vanishes as the temperature 
goes to zero and has a maximum when the temperature 
is of the order of the Thouless temperature T c = E c /kB- 
At elevated temperatures, T » T c , the parameter con- 
trolling the decrease in amplitude of the conductance os- 
cillations is Ec/ksT. This is simply the relative number 
of electrons with energy of order E c . These electrons are 
responsible for the interference effect we arc discussing, 
which is nothing but the "thermal effect" of Refs. [| and 

! 

Now we turn to structures with beam splitters; below 
we show that beam-splitting scattering between differ- 
ent trajectories at the junctions between the mesoscopic 
sample and the reservoirs qualitatively changes the in- 
terference pattern. In this case quasiparticles with low 
excitation energies, E < E c , may contribute — in some 
cases resonantly — to the interference effects causing the 
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conductance to oscillate as a function the superconduct- 
ing phase difference. 



2. Giant conductance oscillations in Andreev interferometers 
containing beam splitters 

Scattering due to potential barriers or geometrical fea- 
tures at junctions between the mesoscopic region and the 
reservoirs qualitativly change the nature of quasiparticle 
trajectories. In particular, a particle reflected from an 
N-S boundary does not necessarily leave the sample for 
the reservoir directly. Instead, it may be reflected by the 
junction and re-enter the mesoscopic region. There is a 
certain probability that such reflections creates low en- 
ergy trajectories that connect the reservoir(s) with both 
superconductors. An example of such a trajectory is 
shown in Fig. |l|. 

The role of beam splitters in Andreev interferometers 
was first payed attention to by Nakano and Takayanagi.u 
A number of other interference phenomena also involv- 
ing quasi-particles at the Fermi energy (zero temperature 
phenomena) has been discussed in the literature. For in- 
stance, Wees et al£3 showed that elastic scatterers gen- 
erate muliple reflections at the N-S boundary resulting 
in an enhancement of the conductance above its classical 
value. In ballistic structures resonant tunneling through 
Andreev energy levels coinciding with the Fermi level was 
predicted in Refs. |l^j2^. For diffusive structures contain- 
ing beam splitters a significant increase of the Aharonov- 
Bohm oscillations of the conductance was shown to ex- 
ist in Refs. ||,||,f9[^|@ Beenakker et alM showed 
that the angular distribution of quasi-particles Andreev 
reflected by a disordered normal-metal - superconductor 
junction has a narrow peak centerd around the angle of 
incidence. The peak is higher than the coherent backscat- 
tering peak in the normal state by a large factor G/Gq (G 
is the conductance of the junction and Go = 2e 2 /h). The 
authors identified the enhanced backscattering as the ori- 
gin of the increase of the oscillation amplitude predicted 
in Refs. p^Jl8|. As a final example we note that it was 
shown in Ref. [21| that the beam splitter violates the "sum 
rule" according to which the conductance in the absence 
of junction scattering is equal to the number of trans- 
verse modes and does not depend on the superconducting 
phase difference. 

All the mentioned interference phenomena involving 
quasiparticles at the Fermi level (E — 0) have the same 
nature for both ballistic and diffusive structures. This 
follows from the complete phase conjugation of electron- 
and hole excitations at the Fermi energy. At the Fermi 
energy even a random- walk-type of diffusive electron tra- 
jectory caused by impurity scattering is completely re- 
versed by the Andreev reflected hole and there is a com- 
plete compensation of phase. In particular the giant os- 
cillations of conductance with phase difference <f> is insen- 
sitive to impurities, as there is a finite scattering volume 



in which the phase gains along the electron-hole trajec- 
tories are completely compensated. 

When the transparency associated with junction scat- 
tering has intermediate values both the thermal ef- 
fect and the resonant oscillation effect contribute si- 
multaneously provided the temperature is close to 
the Thouless temperature. In experiments measuring 
the conductance., oscillations for structures with beam 
splitter JlliHBBErMd the temperature was of the or- 
der of the Thouless temperature or higher, and hence 
both effects could contribute. The effects can be distin- 
guished by lowering the temperature below the Thouless 
temperature, as then the amplitude of the conductance 
oscillations decrease in the case of the thermal effect (it 
goes to zero as the temperature goes to zero) while the 
resonant amplitude of the conductance increases and is 
maximal at zero temperatjujee. Experimental evidence is 
just beginning to appearEjED 

While the role of the termal mechanism has been in- 
vestigated in detail in Refs. [||], f° r the giant conduc- 
tance oscillations the role of intensity of scattering for 
all types of scattering mentioned above (normal [non- 
Andreev] reflection at the N-S boundaries, junction- and 
impurity scattering) remains without a quantative de- 
scription. The objective of this paper is to fill this gap. 

The paper is organized as follows: in Section [□] we de- 
scribe how Andreev interferometers are modelled in this 



work; in Section [II, we develop a resonant perturbation 
theory to find the conductance in the case of ballistic 
transport inside the sample and weak coupling of the 
sample to the reservoirs. In comparison with Ref. [l^ we 
here allow scattering between different conduction chan- 
nels at the two junctions between sample and leads to 
reservoirs. In Section [Fv| we in addition take into account 
the normal reflection that accompanies the Andreev re- 
flection of an electron (hole) at a real normal conductor- 
superconductor interface, and get an explicit analytical 
expression for the conductance of the system as a func- 
tion of the number of transverse channels. For cases when 
it is inconvenient to get analytical results, such as when 
the coupling is not weak, we present some results of nu- 
merical calculations in Section Then, in Section VI 
we relax the condition of the sample being in the ballis- 
tic transport regime and calculate the giant conductance 
oscillations for a diffusive hybrid S-N-S structure using 
the Feynman path integral approach for the trans ition 



probability amplitude. In the conclusions, Section VII 



we discuss the range of parameters for which the conduc- 
tance oscillations can be giant in real experiments. 



II. MODEL FOR AN ANDREEV 
INTERFEROMETER 

In this Section we describe our model for an Andreev 
interferometer. As schematically shown in Fig. ref system 
it consists of a superconductor-normal (semiconductors- 
superconductor (S-N-S) sample coupled to two normal 
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electron reservoirs between which..—, iioltage bias is ap- 
plied. Appealing to experimentsoEacil we neglect scat- 
tering of electrons by impurities inside the sample for the 



time being and return to this point in Section VI. Nev- 
ertheless, the junctions between the S-N-S sample and 
the normal leads to the electron reservoirs inevitably are 
sources of scattering. So, whereas we consider electron 
transport to be adiabatic inside the sample — the current 
being carried in N± channels (modes) — electrons can be 
scattered between different conduction channels at these 
junctions. Taking this into account amounts to a first 
generalization of our earlier treatment of this problem. ta 
In our model the coupling between the sample and the 
reservoirs is controlled by potential barriers (beam split- 
ters, see above) appearing at the junctions between the 
leads from the reservoirs and the sample. We assume that 
in the case of low barrier transparency the approximation 
of a nearly isolated sample is adequate and that channel 
mixing is absent; we shall then study what happens when 



the coupling increases in Section III 




a b a b 



FIG. 4. Schematic picture of an Andreev interferometer 
of the same type as shown in Fig. 1. The full [dashed] arrows 
indicate electrons (e) [holes (h)] moving in the ballistic seg- 
ments 1-5 of the sample. In the model calculation described 
in the text Andreev- and/or normal scattering may occur at 
the two superconducting mirrors (S) and scattering between 
different segments and channels (modes) may occur at the 
junctions marked A and B. 

Another fact ignored in our earlier wor is the pos- 
sible "normal" reflection of quasiparticles by a Schottky 
barrier at the S-N boundaries, a mechanism that would 
compete with the Andreev reflection. When normal re- 
flection is possible the degeneracy of the quasiparticle 
energy levels (Andreev levels) which occurs at the Fermi 
level for certain values of the phase difference between 
the two superconductors is lifted and the "giant" conduc- 
tance oscillations as a function of this phase difference 
is greatly reduced. Despite the experimentaL.faci_.that 
the probability for normal reflection is smalloBeil the 
criteria for how small the normal reflection probability 
must be for the giant oscillations to survive is obviously 
an important question, which we consider in Section [Tv| . 
Below we formulate our transport problem for a general 
case which includes both the possibility for scattering be- 
tween conduction channels at the sample-lead junctions 



and normal reflections at the N-S interfaces. 

It is convenient to divide our model Andreev inter- 
ferometer into five different segments so that the elec- 
tron transport to a good approximation is adiabatic in 
each segment. We then use a phenomenological method 
for describing the two manifestly non-adiabatic junc- 
tion regions [marked A and B in Fig. ref system]. The 
quasi-particle wave functions in the adiabatic segments 
1-5 shown in Fig. [| can be found with the help of the 
Bogoliubov-de-Gennes equation. As channel mixing is 
absent in the adiabatic segments the electron- and hole 
like components of the wave function in the n:th trans- 
verse mode in segment a is 

N± 

u a {x,y) = £ (a&e< e, * + 6« e-*£ )s ) S ink ± (n)y 
n=l 

Nx 

v a (x,y) = g (*£to ik<h)x + €lef h)x ) sink ± (n)y (1) 



Here a^n and ba, n [aa}n and b£j„] are the probabil- 
ity amplitudes for free motion of electrons [holes] for- 
ward and backward, respectively, in channel n and seg- 
ment a of the sample; k±(n) — mt/d,n = 0,1,2,... 
is the quantized transverse wavevector assuming a hard 
wall confining potential, d is the width of the sample, 



[kp - k\(n) ± 2mE/h 2 } 1 / 2 is the electron (hole) 



longitudinal momentum, kp is the Fermi wavevector, E 
is the electron energy measured from the Fermi energy, 
while x and y are longitudinal and transverse coordinates 
in the sample, respectively. Non-adiabatic scattering of 
electrons in the junction regions, see Fig. [|, is described 
by a unitary scattering matrix S connecting the wave 
functions in the surrounding sample segments. Scatter- 
ing at these junctions mix the transverse modes (chan- 
nels) in the adiabatic segments (which here and below, 
for the sake of simplicity, are considered to have the same 
number of open transverse channels). Hence, the scat- 



tering matrix connects c 



(in) 



, (out) 

and Cq , 



which are N± 



-component vectors whose coefficients de- 
scribe the incoming and outgoing adiabatic wave func- 
tions [see Fig. || and Eq. (Q)], 



■.{out) 



(2) 



We assume the coupling matrix S to be symmetric with 
respect to the left and right sample segments [labeled 2 
and 4 in Fig. Therefore the matrix can be written as 



Sll S12 S12 
§12 §22 §23 
§12 §23 §22 



(3) 



where s a p are N±_ x Nj_ matrices which mix the con- 
duction channels when an electron (or hole) is trans- 
ferred from the (3 - to the a segment. Electrons and 



G 



holes are, however, not mixed. The elements of 
(n,m — 1,2, ...N±) are the probability amplitudes for 
an electron (or hole) in the m-th channel of the /3-section 
to be transferred to the n-th channel of the a-section. 
We assume that scattering of an incident quasiparticle 
at the junction causes transmission of the electron into 
each of the N± open transverse channels with a prob- 
ability, which is of the same order of magnitude for all 
channels. This implies that the matrix elements of the 
matrices §12, S22, S23 — 1 and su — 1 are of order 1/yNx- 



(in) 
C 2 



(in) 
C3 



,(out) 



(out) 



,(out) 



(in) 



1 



FIG. 5. Detail of the junction (A, B in Fig. 4) coupling 
the reservoirs via leads to the normal part of the system. A 
scattering matrix connects the amplitudes of incoming- and 
outgoing quasiparticles, see text. 

We choose to parametrize the S matrix in a way such that 
there is no channel mixing if the sample is completely de- 
coupled from the reservoirs. This coupling is determined 
by the elements of the matrix §12 which are the prob- 
ability amplitudes for electron (hole) transitions from a 
lead (segments 1 or 5) to the sample (segments 2, 3 or 
3,4). In order to describe the strength of the coupling we 
introduce the parameter e r and write 



S12 



C-r 



1/2 



S12 • 



(4) 



The scattering matrix (||) has to be unitary — a require- 
ment that leads to five relations between the eight matri- 
ces s aj g, see Eqs. ( f79f|83| ) in Appendix A. (Note that each 
matrix has an independent Hermitian and anti-Hermitian 
part). We are thus left with three undetermined matri- 
ces, which we choose to be §12 and the anti-Hermitian 
part of §22 (the last choice is made only for the sake of 
calculational convenience, see Appendix A. We assume 
the Hermitian and anti-Hermitian parts of §22 to be of 
the same order in the parameter e r since in the general 
case they are connected by a Kramers-Kronig relation. 
It follows from the unitarity conditions that the matrix 
elements of §12 are of order 1 / . Hence the matrix 
elements of S12 are of order unity. 

The conductance of our model system is in the limit of 
vanishing bias voltage determined by the Landauer for- 
mula as modified by Lambert for a system with Andreev 



reflection: 



J: 



2(R A R' A — T A T' A ) 



R A + R A 



(5) 



Here e is the electron charge, h is Planck's constant and 



N ± N± 

TA = J2ri A \R A = ^ P ^ 

k=l k=l 

N ± N ± 

to = e^ 0) ^ = E^ 0) - 
fe=i fc=i 



(6) 



] is the probability for an electron 



In Eq. (f§ 

approaching the sample in the fc-th transverse channel of 
the left lead to be transmitted as an electron (hole) into 
any of the outgoing channels of the right lead. The quan- 
tity p^ [p^] is the probability for the same electron to 
be reflected as an electron (hole) into any outgoing chan- 
nel in the same left lead it came from. Similarily, t 



'(<>) 



[t^ '] and p' fc w [p'^'} are normal (Andreev) probabili- 
ties for an incoming electron from the fc-th transverse 
channel of the right lead to be transmitted as an electron 
(hole) into any outgoing channel of the left lead and to 
be reflected as an electron (hole) back into any outgoing 
channel of the right lead, respectively. 

In order to proceed we have to solve the matching 
equations for the adiabatic wave functions in sample and 
leads. The matching problem under consideration is il- 
lustrated in Fig. U where solid and dashed arrows symbol- 
ically show electron- and hole plane waves moving to the 
right and to the left, respectively. The coefficients a^ 6 
and ha are 7V_L-component vectors, the components 
of which are the probability amplitudes aa'n^ and ba,n , 
see Eq. (Q). Matching the wave functions at the junctions 
using Eq. (||) one gets the following set of equations for 
these amplitudes: 



M) _ 5 _(e,fo) 

L>1 — 



S11&5 



(e,h) 



te.h 

si 2 (7, a. 



§12 U c 



->(e,fr)i (e,h) 



(7) 



and 



u( e >'0 a a ( e > h ) 
d 2 - s 22 a 2 

(e.h) _ £ (e,h) 



(e,h) 



a 



S23&2 



u 



(e,h)i(e,h) 



§22^; 
§23^; 



S22bg = s i2 a^ 

(e,h) (e,h) 
3 A 3 
(e,h) (e,h) 
a. 



■ §23 V t 



{e,h)\,{e,h) 
i D 4 



= Sl2& 5 

(e, 

(8) 



(« 



(e,h) 



Here the diagonal matrices U a simply keep track of 
the phase gained by electrons and holes during their free 
motion across segment a. The diagonal matrix elements 
are 
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ui e) (n) = u ( £\n) = e - ik « )l °, ri, = 1, 2, . . . N x 

(9) 

where l a is the length of section a in Fig. || 

The set of equations (Q) and (||) must be supplemented 
with boundary conditions at the N-S interfaces. In the 
general case when both Andreev- and normal reflections 
at the N-S boundaries are possible the boundary condi- 
tions are 



electrons perturbed by a coupling of the sample to the 
reservoirs are expected to be found (Hvf/L is the char- 
acteristic spacing of mode energy levels near the Fermi 
energy). 

By solving the matchi ng eq uations (|J), (| io[), and (|1] 
for e r = 0, using Eqs. ( |l08| ), ( |TTo| ), and M of Ap- 
pendix B one recovers the following expressions for the 
low energy spectrum of electron-hole quasiparticles in the 
normally conducting semiconductor sandwiched between 
two superconducting mirrors 



W e u 9 r, T I A C U 9 „ 



JV 



Jh) _ [ (1)* iSk(n)h h (e) , (1)* p -i2fc< h >(n)i 2 ,(h)l / lf) \ 

"2,n — e ' A C U 2,n ' N C U 2,n \ lyJ ) 

for the left (first) boundary [Sk = k^ — k^} and 



°4,n — e 



(2)„(e) 
"AT a 4,n 



J A,n 



J^2 



'A u 4,n 



4,n 



JV 



(11) 



for the right (second) boundary, see Fig. ^. The probabil- 
ity amplitudes for normal- and Andreev reflections at the 
N-S boundary are given bye 4 *r/v and e 1 *^. It follows 

For convenience explicit ex- 



that 



„U.2)|2_ L |„.(l,2)|2 
JV 



by ( 
1§ 



pressions for these quantities in terms of the complex or- 
der parameter of the superconductor and the reflection- 
and transmission probability amplitudes of the normal 
barrier at the N-S interface are given in Appendix B. 

The equations (@) , (ft) , (0) , and (|ll|) together with the 
conductance formula (pt) form a complete set of equations 
that permits us to find the conductance of the system 
under consideration. 

In the next section we discuss the non-adiabatic scat- 
tering of electrons at the junctions and present analytical 
formulae for the case of a weak coupling of the sample 
to the reservoirs (e r ^ 1) and numerical results of com- 
puter simulations in the general case. The role of the 
non-Andreev (normal) reflection at the N-S boundaries 
is discussed in section IV. 



III. ROLE OF SCATTERING AND MODE 
MIXING AT THE POINTS OF COUPLING TO 
THE RESERVOIRS 

We start our analysis by assuming a weak coupling be- 
tween sample and (leads to) reservoirs. In this case the 
parameter e r introduced in Eq. (Q) is much smaller than 
one. It is convenient to develop a qualitative understand- 
ing starting from the so called Andreev levels that form in 
the isolated sample when e r is strictly zero. We consider 
values of <fi near odd multiples of 7r for which Andreev 
levels will appear at the Fermi energy (0 = (f>2 — (pi; 4>\ 
and 4>2 are the phases of the gap functions in the left 
and right superconductor, see Fig. |]). We concentrate 
on energies in a narrow interval AE ~ e r hvF/L around 
the Fermi energy, whithin which the quantum states of 



tt(2Z + l)±<f>= (jk F - fex(n) 2 + 2mE/Ti 2 



(12) 



kp — kj_(n) 2 — 2mE/h 2 ) L, / = 0,±1,±2 ) 



Here m is the electron mass, L = 1-2 + h + h is the length 
of the normal part of the sample, and all reflections are 
assumed to be of the Andreev type. The phase <fi comes 
with a plus- or a minus sign in Eq. ([l2|) depending on 
whether the electron-hole excitations move as electrons or 
holes when going from the left to the right S-N interface. 

Expression (|12|) for the spectrum tells us that when 
4> = n(2l' + 1) there is one Andreev state at the Fermi 
level for each transverse mode (index n) simultaneously, 
i.e., the energy of the state whose quantum number I 
associated with the longitudinal motion equals I' coin- 
cides with the Fermi energy irrespective of mode number 
n. Therefore, the degeneracy of the energy level at the 
Fermi energy (E — 0) is given by the number of open 
transverse modes N±, whenever (f> equals an odd multi- 
ple of 7r. This results in a giant probability for resonant 
transmission of electrons from one reservoir to the other. 
The amplitude of the carresponding conductance oscil- 
lations, AG oc iVie 2 /ML3, is therefore much larger than 
the conductance quantum. 

A finite coupling of the sample to the reservoirs (which 
is of course necessary for a current to be observable) si- 
multaneously results in a broadening and a shift of the 
Andreev energy levels. The former effect is due to quasi- 
particle tunneling from the sample to the reservoirs after 
a finite time, the latter is due to mixing of the transverse 
modes that inevitably accompanies a finite coupling. Be- 
low we show the broadening and the shift to be of the 
same order in the transparency of the barrier connect- 
ing sample and reservoirs. The result is a broadening of 
the peaks of resonant sample conductance but their giant 
amplitude remains. This is beacuse a Breit-Wigner type 
of resonance is broadened without loss of amplitude when 
the coupling is increased. It turns out that the broaden- 
ing of each state tends to compensate the shifting around 
of the energies of previously degenerate states. Read- 
ers who are not interested in technical details may want 
to turn directly to Eq. (p9|), which expresses this result. 
Results of numerical calculations presented in Section |v| 
show this picture to hold up to a value for e r which is 
about half its maximum value. A further increase of the 
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coupling results in a large decrease of the amplitude and 
increase of the broadening of the peaks. 

In the weak coupling case the set of equations (fjj), (||), 
(|Io|), and (|ll| ) which determines the transmission prob- 
ability amplitudes can be solved analytically by pertur- 
bation theory in the small parameter e r . Below this per- 
turbation theory will be developed. 

As shown in Appendix A, all the matrices which de- 
scribe scattering of electrons and holes inside the sample 
(J22 and S23 — 1) are proportional to e r if the coupling 
matrix s 12 is proportional to ^/e^ and if e r -C 1 . Hence 
it follows that the set of equations (0-^l|) can be written 
in the form 



(W(E) 



e r n) |H) = 



1/2 



IK) 



(13) 



The vector |H) has 12 N± unknown components, viz. 
af' h \ b^' l) for i =2, 3, and 4. The vector |K) has 
AN± known elements, a.^ for k =1 and 5, and 8N± 
elements which are zero. The matrix ~W(E) has 12A^_l 
block matrices along the diagonal with non-zero elements 



>rmW n n(E) 



(14) 



where a, f3 — 1, 2, ... 12; n,m — 1 . . . N±. The matrix 
wW(.E) has been obtained for the n-th fixed transverse 
mode in the absence of coupling (e r = 0) by matching the 
electron- and hole components of the wave functions at 
the N-S boundaries using Eqs. @ and (|l|) and at the 
junctions coupling the sample to the electron reservoirs 
using Eq. (||) for fixed channel number n and e r = 0. 
The matrix f2 has elements describing mixing between 
modes. The explicit forms of operators W, fl and vec- 
tors |H), |K) are straightforwardly found by comparing 
Eqs. (|) and ' 

In order to use the resonant perturbation theory we 
have to consider some properties of the unperturbed sys- 
tem relevant to our problem. It is straightforward to see 
from (|lj) that the determinant of the matrix W(-E) can 
be written as a product of N± factors 

DetW(E) = Detw (1) (E) x . . . x Detw^ (E) (15) 

and that its value iz zero at any eigenvalue E = E n ; 

(n)\ 



the energy has to be much smaller than the distance be- 
tween neighbouring energy levels corresponding to quan- 
tization of the longitudinal motion of electrons, that is 



e r hvF/L <C Tivf/L 



(18) 



Here we develop the perturbation theory for a general 
case in order to use the results also in the next section. 
Therefore in order to find the correct zero order wave- 
function the vector |H) must be taken as a superposition 
of the Nr states inside the resonant region (Nr is possi- 
bly but not necessarily smaller than N±) 



IH) 



£ 

n=l 



7 „|# ) > + |H 1 ). 



(19) 



The summation in (|T^) goes over the Nr transverse 
modes inside the resonant region, which extends over an 
interval of order e r hvp/L on either side of the Fermi en- 
ergy; |Hi) is a small addition oc e r . The unknown coeffi- 
cients 7„ should be found with the solvability condition 
of the equation for |Hi) that is readily available from 
Eq. (|13|) m the linear approximation in e r <l: 



Nr 



W(£)|Hx) = -J2ln[W'(E n j)(E - E n j) - Erfl}]^ 



(«)\ 



(20) 



Here the superscript "prime" indicates derivation with 
respect to energy E. When obtaining Eq. ( p0| ) we used 
the inequality (|l^) and expanded ~W(E) in a Taylor series 
around E n ^i (with the restriction \E — E n> i\ <C hvp/L) in 
every term of the sum and took into account Eq. (Eq). 

Multiplying both sides of Eq. ( po|) from the left by bra- 
vectors | [which can be determined from the equa- 

tion (4 m) \W(E rn ) = 0] one readily gets the solvability 
conditions for Eq. ( p0| ) that determines the coefficients 
7„. In this way we obtain the main equation which has 
to be solved in order to get 7 n ; these coefficients, accord- 
ing to Eqs. (Q) and ([l9|), determine the probability of the 
resonant transmission of an electron from one reservoir 
to the other via the sample) : 



N n 



of the unperturbed system. The eigenfunctions \ip\ ) of \i(E — E )W' S — e fl 1 

the unperturbed problem satisfy the following equation ^ n ' n mn e ' mn * ' 



n=l 



1/2 



K„ 



W(E n<l )\^ n) ) = 



(16) 



Developing the perturbation theory we assume the fol- 
lowing inequality to be satisfied. 



Here we have used the short hand notation 



X F < d < L , 



(17) 



where Xp is the de-Broglie wave length (Fermi wave 
length) of the electron, while d and L is the width and 
length of the sample. We note that the perturbation of 



W' n - -i{ipW\W\E n )\i>W) 

fi ro „ = <V (m) |%W) 



(21) 

(22) 
(23) 
(24) 
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We have also dropped the subscript I as we have as- 
sumed it does not change under the perturbation con- 
sidered. Using Eq. ( |l6| ) for e r = it is straightforward 
to calculate W n and show it to be real, i.e., the Hermi- 
tian and anti-Hcrmitian parts of the coupling matrix 
provide broadening and shift of the energy levels of the 
sample, respectively. In our analysis of Eq. (pl| ) we [Con- 
sider the matrix elements f2 nm to be of order unity.O It 
is then easy to see that far from the resonance, where 
Tivf/L ^> \E — E n \ ^> e r TiVF/L, the first term on the 
right hand side of Eq. (21) dominates and one gets 



In 



1/2 



iW' n (E-E n ) 



(25) 



Knowing j n we may calculate |H), which contains the 
coefficients and 64, from Eq. ([lj]). By using Eq. (Q) 
the probability of transmission of an electron (hole) from 
one reservoir to the other is 



In the range of resonant energies, \E 
the amplitudes 7, 



therefore, the transmission probability amplitudes 



(26) 

E n \ < e r hv F /L, 
are much larger and, 



5,m 



'12 



*3,n 



(27) 



obtained from Eqs. (fij), ([To|) , and (||) are independent of 
e r . Note that „„a fc (o)(e,h) _ ^ 1 — 1; 



l 3,n 



and 



are the known ampli- 



tudes of the wave function of the electron (hole) in the 
n-th transverse mode in sample segments 3 and 4 when 
isolated from the reservoirs. Hence c4°n 6 and bf^'^ 
are of order unity and it follows that the probability for 
an electron in the m-th transverse mode of segment 1 
the lead from the left reservoir — to be transmitted to 
any of the N± transverse modes in segment 5 — the lead 
to the right reservoir, (see Fig. ||) — via the sample is 



rp(m ) _ \ " 
1 (e,h) - 2^ 



j,(e,h)|2 
J 5,m I 



EE 

m—1 n—1 



mn\2 I (rrio) |2 



'12 



17, 



Nr 
N ± 



(28) 



The last similarity relation follows since J\=i I s i2 1 
e r and since in the resonance region, according to 

Eq. @, |7i mo) | - l/VN±e^ (to see this note that 
Eq. ( pq ) is valid up to the resonant region where \E — 
E n \ ~ e r hvF/L; the superscript mo indicates that the in- 
coming electron in segment 1 moves in mode number mo). 
Therefore in accordance with the Landauer-Lambert for- 
mula (^), the order-of-magnitude conductance in the res- 
onant region of a system with N± incoming electrons is 



N± 

1 (e,h) 

m = l 



e 2 



(29) 



while off the resonance G ~ efNne 2 /h [cf. Eq. (26)]. 

Since at zero temperature the energy of the incom- 
ing electrons coincides with the Fermi energy, resonant 
transmission occurs in the vicinity of <f> = ir(2l + 1), I = 
0, ±1,±2, the width of the resonance being of order 
e r <C 1. If reflections from the N-S boundaries are only 
of the Andreev type it follows that Nr in (^9|) is equal 
to R±. In this case the conductance oscillates with <fi, 
the amplitude of the oscillations being proportional to 
the total number of the transverse modes N±. In the 
above analysis, for the sake of simplicity, we assumed the 
number of transverse modes inside the sample and the 
leads to be equal but it can easily be shown that if these 
numbers are different the conductance is proportional to 
the smallest one. 

As demonstrated in this Section, for the many-channel 
case with mixing of transverse modes at the junctions our 
analytical approach permits us to estimate the conduc- 
tance in the region far from the resonance. It is also pos- 
sible to find the width of the resonant peak and its height 
(i.e., the amplitude of the conductance oscillations) but 
it does not permit us to find the fine structure of the 
resonant peak as it is determined by the set of A^_l 3> 1 
algebraic equations of Eq. ( pl| ) . Here we consider instead 
the fine structue of the resonant peak using the most sim- 
ple model of a one-channel sample weakly coupled to the 
reservoirs. In this case calculations of the conductance 
in the vicinity of the resonance (S<p = \<f> — tt| -C 1) give 
the result 



G = 



2e 2 



(4 7 ) 2 



h [(2 7 )2 + ((ty)2 



[(27 cos kl 3 



>Hinfcf3) a + (ty) 2 ] 
(30) 



(k is the electron wave number, 
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is the distance be- 
tween the junctions, 7 = |si2| 2 ~ e r "C 1). It follows 
that there is a dip in the middle of the resonant peak 
(which appears due to an interference between the wave 
functions of the clockwise and counter clockwise motions 
of the quasiparticles) . When 5(f> = the conductance is 



2e 2 2 
G = — — cos M3 
h 



(31) 



and hence is it goes to zero for certain values of the wave 
number k\ the resonant peak is split into two peaks. 

In the many channel case every mode has its own lon- 
gitudinal momentum, and the conductance being a sum 
over the channels is self-averaged with respect to momen- 
tum. Such an averaging of the conductance in Eq. ( |30| ) 
followed by a multiplication by the number of transverse 
modes gives as a result for the conductance, 



2e 2 



" V - 2 " K27) 2 



(2 7 ) 2 + 3(^) 2 



(32) 



This result tells us that there is a dip in the middle of the 
resonant peak with a depth of 1 /9 of the height of the 
resonant peak. Equation (p3) is valid in the absence of 
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transverse mode mixing. Numerical calculations of the 
conductance for the general case of the transverse mode 
mixing also show such a dip in the middle of the resonant 
peak (see below). 



IV. INFLUENCE OF NORMAL QUASIPARTICLE 
REFLECTION AT THE N-S BOUNDARIES ON 
THE GIANT CONDUCTANCE OSCILLATIONS. 

In experiments a typical N-S boundary is an interface 
of two different conductors, resulting in two-channel re- 
flection of electrons at the N-S boundary that is an in- 
cident electron is reflected back remaining in the state 
of an electron- like excitation with probability |r/v| 2 (the 
normal channel) and in a state of a hole-like excitation 
with probability \ta\ 2 = 1 — Vn\ 2 (the Andreev channel). 
In the general case of nonequivalent normal barriers at 
the NS boundaries the quantized energy levels of an S- 
N-S system are repelled from the Fermi level and the de-i 
generacy is lifted. However, we know from experimentsEa 
that a situation with a low probability for non- Andreev 
(normal) reflection can be realized in practice. Therefore 
it is important to derive a criterion for how low this prob- 
ablity for normal reflection must be to preserve the giant 
conductance oscillations. In this Section we discuss the 
role of the normal reflections for the oscillations of the 
conductance in a ballistic S-N-S system with combined 
Andreev and normal reflections at the S-N boundaries. 



A. Normal reflection from two identical barriers at 
the N-S interfaces 

We start from the case of a sample isolated from 
the reservoirs, the geometry of which is presented in 
Fig. ||, and assume the reflection properties at the two 
NS boundaries to be identical. Matching the wave func- 
tions of the electron- and hole-like excitations at the N-S 
boundaries using Eqs. ( |l0| ) and (|ll]) gives as a result the 
following spectral function, 



cos<^>_ — |rjv| cos tp+ + \ta\ cos< 



(33) 



Here <p_ = 2mEL/h k n , ip+ = 2k n L, the parallel compo- 
nent of the wavevector is k n = [k 2 F — fcj_(n) 2 ]^l/2) where 
k±(n) is the projection of the wavevector on the N-S 
boundaries, n labels transverse modes and cj) is the phase 
difference between the order parameters in the two su- 
perconductors. 

For energies E small compared to the energy gaps in 
the superconductors the equation Q n — determines the 
discrete Andreev energy levels of the system. This rela- 
tion can be rewritten as 



E nd = [tt(2Z + 1)± 

arccos (\ta\ 2 cos0 — \tm\ 2 cos<^ + )] 



(34) 



where the longitudinal and transverse quantum number 
is I = 0, ±1, ±2, ... and n = 1, 2, N±, respectively. 

In the absence of normal reflection at the N-S bound- 
aries (rjv =0) Eq. (34) reduces to Eq. (|l2|) and the energy 
level at the Fermi energy is iVj_-fold degenerate for values 
of </) that corresponds to odd multiples of tt. This is the 
case described in the previous Section. For the symmet- 
ric case of equivalent barriers at the two N-S boundaries 
the normal reflection lifts this degeneracy at the Fermi 
energy, as can be deduced from Eq. (|34|). We show below, 
however, that the lifting of the degeneracy is restricted 
in the sense that the amplitude of the giant conductance 
oscillations remains proportional to N± . 

We begin with a qualitative argument and neglect as a 
first step the quantization of the transverse momentum. 
Hence we consider k±{n) to be a continuous variable 
(k±(n) — > k±). Within this approximation the spectrum 
Ei(k±) and the wave functions \l,k± > of a quasiparti- 
cle are characterized by one discrete quantum number I 
associated with the longitudinal quantization and by one 
continuous variable, the transverse wave vector k±. As 
can be seen from Eq. ([34|) energy levels are at the Fermi 
energy [E = 0) if two conditions are satisfied, viz. 



and 



7r(2s + 1) 



lf + = 2kL = 27TQ0 



(35) 



(36) 



2mL 



where k (we have dropped the subscript n) now is a con- 
tinuous variable; sq and go are integer numbers. It follows 
that in the absence of transverse momentum quantization 
the symmetric barriers at the N-S boundaries do not com- 
pletely remove the degeneracy of the energy level at the 
Fermi energy. The extent of the degeneracy depends on 
the number of transverse wavevectors [cf. Eq. ([56])] for 
which the equation ( |34| ) is satisfied. This number is de- 
termined by the largest possible value of qo, which will 

be estimated below. 

From its definition one notes that k = \/k'j? — fc^" 
varies between zero and Uf and hence from ( |36| ) one 
concludes that < qo < kpL/ir. This implies that 
the maximum value of qo, let's call it Nq, is of order 
kpL ^> 1. Therefore, whenever cj) — ir(2s + 1), there 
is a degenerate energy level at the Fermi level with de- 
generacy ~ Nq. The number of states through which an 
electron can be resonantly transmitted from one reser- 
voir to the other is even greater, however. This is be- 
cause the width of the energy levels broadened due to 
the coupling of the sample to the electron reservoirs is 
5E ~ e r hvp/L and all the quantum states inside this 
range of energy resonantly transfer reservoir electrons 
through the sample. In order to determine the num- 
ber of states within this energy range we estimate the 
total width of the intervals in the fc-space around points 
k = irqo/L inside which wave functions \l,kx > of the 
system correspond to energy levels inside this range of 
energy, Ei(k±) < e r hvp/L. We do so by expanding the 
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cosines in Eq. ( p4| ) in a Taylor series in the small de- 
viations 5k and 5E = e r hvp/L near one of the points 
where the cosines are equal to unity [these points are 
determined by Eqs. (|35| ) and pf])]. Employing the sum 
rule \rA\ 2 + |rjv| 2 = 1 one can show that the energy lev- 
els are inside the resonant range E < 8E — e r hvp/L if 
5k < e r /\rN\L, assuming |rjv| 3> &r- Multiplication by 
the number Nq of such intervals gives the total range of 
the "resonant" momenta as 



Afc 



£r < \tn\ 



(37) 



A similar analysis shows that if e r > |rjv | all N± trans- 
verse modes take part in the resonant transition; the os- 



cillations disappear if \ta\ 



211/2 



< £r 



Now we go one step further and take the trans- 
verse quantization into account. In the limit 1/JVj. <C 
e r -C 1 the quantized values of momentum k n — [k 2 F — 
(mr/d) 2 ](l/2) are almost evenly distributed between zero 
and hp. Hence it follows that the probability for a 
transverse mode to be inside the resonant interval Afc 
is P = Afc/fcf = e r /|rjv|- Therefore the total num- 
ber of transverse modes inside the resonant region Afc is 
N R ~ N±P = N ± e r /\r N \. From here and from Eq. @ 
it follows that the maximum conductance (when elec- 
trons are resonantly transmitted through the sample) is 



Grain ^ 
Gmax 2 



M 2 + e 2/4 

M 2 + e 2 /4 



(42) 



Hence it follows that 

el/A 



G„ 



G r , 



2\r A 
1, 



|rjv| < £r 
|rjv| ~ \ta\ 



(43) 



In a situation when |rjv| ~ \ t a\ the amplitude of the 
conductance oscillations is greater by a factor N±e r ^> 1 
than in the absence of the superconducting mirrors. If 
|rj\r | < £r the amplitude of the conductance oscillations 
is oc Nj_. 

In the above analysis we considered the case of equiva- 
lent boundary potentials, so that the probabilities of nor- 
mal reflection are equal at the two N-S interfaces. When 
these probabilities are not equal, the energy levels never 
reach the Fermi energy and resonant transmission occurs 
only if the asymmetry is not too large. Below we analyse 
the situation of non-equivalent N-S boundary potentials. 



B. Normal reflections from non-equivalent N-S 
boundary potentials 



G r , 



oc N 



2e 2 



(38) 



Anal ytic al calculations presented in Appendix C, see 
Eq. (12C), show the conductance of a sample with sym- 
metric N-S boundaries (i.e., boundaries with equal prob- 
abilities of normal reflection) to be 



G~N±- 



2e 2 



^(l + lMPcos^ + e 2 ^) 2 -!^! 4 



(39) 



As is evident from Eq. (p9|), the maximum conductance 
occurs when <fi = n(2l + 1), which is when energy levels 
line up with the Fermi energy and, therefore, resonant 
transition of electrons from one reservoir to the other via 
the sample takes place. Using Eq. ( |39| ) it is straightfor- 
ward to sec that the maximal conductance is 



Gr, 



2e 2 



h V\tn\ 



(40) 



If |r/v| <C e r we have the giant conductance oscillations 
predicted in Ref. [l^. If |r/v| 3> £r the maximal conduc- 
tance is determined by Eq. (|38|). The minimal conduc- 
tance — occurring when <f> — 2nl — when we are off 
resonance, is 



G r , 



Na 



2e 2 



e 2 r /2 



VM 2 +e 2 /4 



(41) 



The ratio between the maximal and the minimal conduc- 
tances is therefore 



Matching of the wave functions of the electron- and the 
hole-like excitations at two non-equivalent N-S bound- 
aries results in a spectral function of the form 

r\ I Mil ( 2 )l , I Mil ( 2 )l J. 

Q n = cos <yj_ - \ry I \ry \ cos ip + + \r\ \ \r\ \ cos <p 

(44) 

and the energy levels of the system are determined by 
solutions to the equation 



cos 2mEL/h 2 k n 



= \r$\\rP\ 



cos 2fc„L 



|r^||r^| cosi 



(45) 



Here r N ' and r A are the probability amplitudes 
for an electron to be normally and Andreev reflected, 
respectively, at the left (1) and right (2) boundaries; 

l r !v 2) ! 2 + \ r A ,2) ? = L As follows from Eq. Q, if r$ 



(2) 

and r N are different there is an energy gap in the spec- 
trum around the Fermi energy since the maximal value 
of the right side of Eq. ([lij) is smaller than unity and 
hence there is no energy level at the Fermi energy for 
any <j). For a weak asymmetry between the boundaries, 
<5r/v = |tjvi — rival 1, the maximal value of the right 
hand side of Eq. (H) differs from unity by an amount 



SM = (5r N f 



(46) 



Hence it follows that resonant transmission of electrons 
occurs only if 5r^ < e r - Analytical calculations carried 
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out for the general case in Appendix C shows the con- 
ductance to be 



2e 2 



yj(l + \$ 1 \ r f | cos cp + £ 2 /2) 2 - | r W p \ r f |: 

(47) 



It follows from Eq. ( |47|) that the maximal conductance 
is, with (5r\4 = - r\ '|, 

2e 2 



V /^ + e 2(l_| r W|| r W|)+ c 4/4 



(48) 



Therefore the giant oscillations are of the same kind as 
described above if Sta < e r , but the maximal value of the 
conductance decreases with increasing Sta', when 8ta 3> 
e r the maximal value of the conductance is 



Gr, 



h Sta 



(49) 



V. NUMERICAL CALCULATIONS 

In the range of parameters where e r and hence the 
coupling between sample and reservoirs is not small the 
approximations used above are not valid and the set of 
equations Eq. (Q) must be solved exactly. In order to find 
the largest value of e r for which the conductance oscilla- 
tions are giant, and to find the dependence of the con- 
ductance on parameters of the system we have resorted 
to numerical methods. We have solved the problem for 
different coupling strengths (from 20% to 100% of the 
largest value of e r for which the scattering matrix S of 
Eq. (^) is still unitary; see below), for a varying number 
of transverse modes N± (from 5 to 40) , and for different 
values of the phase difference <f> between the two super- 
conducting condensates (from to 2tt). 

To calculate the conductance of our system we use the 
Lambert formula. The transmission and reflection am- 
plitudes are calculated by matching the waves. Our task 
is to find the probability amplitudes for b\ and 65 for 
quasiparticles going into the reservoirs as functions of pa- 
rameters of the system and of the amplitudes a± and 05 
of quasiparticles approaching the sample from the reser- 
voirs. One parameter is the number of modes N±_ > 1, 
which we relate to the the width of the normal conductors 
(assuming a 2D system) as 



W = (Nj_+0.5)\ F /2 . 



(50) 



The matching of amplitudes at the left (1) and right (2) 
junctions are performed using the scattering matrix of 
Appendix A 



(51) 







( ax 




-4 




a-3 J 




\ b 3 




Sr- 




(52) 



First we eliminate ai and 64 by expressing them in terms 
of &2 and 04, 



a 2 = U2Riu 2 b2 = aib 2 

64 = U4 : R r U4 : a4 = a r a,4 



(53) 



In the next step we eliminate 04 and 62 and to proceed 
we first define 



and then 



ft = a t (l - S221U1) 1 

$ r = & r (l - S^rCtr)' 1 



in = hn + hsifhhii 

121 = $331 + $32lf3i§23l 
Tlr = $21r + S23rPrS31r 
l2r — $22r + S23rPrS32r 



(54) 



(55) 



Using these quantities we conveniently can find the fol- 
lowing expressions for 03 and 63, 



a-3 = Til ai + l2ib3 

63 = (1 - U372r"372i) _1 "3(7lra 5 + ^rU^llCLi) 

Finally we can calculate 



(56) 



h = (sh; + SX21$IS2U)0-1 + {$131 + Si2l$l$23l)b3 

b5 = («llr + Si3r$rS3lr)a.5 + ($12r + Si 3r $ r S 3 2r)u3a3 

(57) 

The studied system is symmetric in the sense that the 
two scattering matrices connecting the sample and the 
reservoirs are equal and the probability of normal reflec- 
tion is the same for both superconducting mirrors. These 
symmetries makes further simplifications possible. 

According to the discussion in section [II (see Ref. |44| ) 
the scattering matrix S in Eq. (||) can be taken to be a 
random matrix. For our numerical calculations we deter- 
mine it as described in Appendix A. 

The scattering matrix describing coupling and mode 
mixing at the junctions have been realized in two differ- 
ent ways. First by assigning random numbers to its ele- 
ments. Here a critical value e c of the coupling strength 
e r was found in the sense that the scattering matrix was 
non-unitary for e r > e c . We find it convenient to define 
a new parameter e = e r /e c , which can be varied between 
and 1. The results from these calculations are shown 
in Figs. Every point is an average of 10 realizations 
of the random scattering matrix. The spread in conduc- 
tance was 4e 2 //i when normal reflection was absent at 
the N-S interfaces and 2e 2 /h when the normal reflection 
probability was at its highest studied value. The position 
of the peak was not seen to change for different realiza- 
tions. The critical value of the coupling was in this case 
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determined by the highest eigenvalue. This gave as a re- 
sult that only some modes where strongly coupled in the 
limit of high e. 
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FIG. 6. Difference 8G between maxima and minima in the 
conductance as a function of the parameter e, which char- 
acterizes the junction scattering matrices (realized by first 
method mentioned in text). Results are plotted for varying 
probabilities R — rjv | 2 for normal reflection at the super- 
conducting mirrors. The number of transverse modes are 
N± = 40. The results agree with the weak coupling limit 
calculated analytically, see Eq. (58). 



parametrized by the coupling parameter e. The result 
of these calculations are shown in Figs. ^-|o[ The only 
parameter to be changed in order to get different realiza- 
tions of the random scattering matrix was an angle (fa 
which only changed the position of the resonant peak. 
For zero angle the shape of the peak is seen in Fig. [ll]. 
The number of open modes are in this realization equal 
to the size of the matrix as all eigenvalues have an am- 
plitude of unity. 

The main result from the analytical calculations to be 
compared with the numerical results is G max — G m i n . 
This is in general approximately equal to G max - From 
equation (p8l) we get 



G r . 



oc N± 



2e 2 



(58) 



which agrees with numerical results when e r < |nv|- 

The first realization of the random scattering matrix 
has been found to describe the weak coupling case as the 
observed peaks were narrow even for e = 1. The second 
realization with a separate matrix mixing modes gave the 
possibility to study weak and intermediate coupling and 
the amplitude of oscillations were seen to diminish when 
e was increased, see Fig. O. 



15.0 




FIG. 7. Difference SG between maxima and minima in 
the conductance as a function of the probability R — \r^\ 2 
for normal reflection at the superconducting mirrors (same 
data as in Fig. 6). The results agree with the weak coupling 
limit calculated analytically, see Eq. (58). 
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FIG. 8. Difference SG between maxima and minima in 
conductance as a function of coupling e (scattering matrix re- 
alized by second method, see text). The probability for nor- 
mal reflection |rjv| 2 = 0. The results agree with the analytical 
results in the weak and intermediate range of coupling where 
the resonant peak is proportional to the number of channels. 



The second type of realization jaf the scattering matrix 
was done with Biittiker matrices^ describing a coupling 
mode by mode between sample and reservoirs. In this 
case an additional unitary matrix was used, which only 
mixed the modes, see Appendix A. Both matrices were 
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FIG. 9. Difference SG between maxima and minima in 
conductance as a function of coupling e. The number of 
open transverse modes are N± = 30. The results agree with 
Eq. (58) for weak coupling. For strong coupling the giant ef- 
fect vanishes in a 2D sample according to the discussion in 
Section I. 
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FIG. 11. The resonance peaks at zero probability for nor- 
mal reflection at the N-S boundaries. The results are from nu- 
merical calculations with N± = 30 and for different strength 
of coupling e = [0.2,0.4,0.6,0.8,1.0] where the most narrow 
peak is for weakest coupling t = 0.2. Note that the amplitude 
of oscillation is much larger than the conductance quantum 
even for e = 1. This is because quasiparticle waves may pass 
the junction in our realization of the scattering matrix even 
if e = 1. 




0.0 0.2 0.4 0.6 0.8 



FIG. 10. Difference SG between maxima and minima in 
conductance as a function of normal reflection probability 
|rjv| 2 (same data as in Fig. 9). The number of open trans- 
verse modes is N± — 30, results for different strength of the 
coupling are shown. The results agree with analytical calcu- 
lations. The |rjv| 2 -dependence agress with Eq. (58). 



VI. GIANT CONDUCTANCE OSCILLATIONS 

FOR A DIFFUSIVE NORMAL SAMPLE - 
FEYNMAN PATH INTEGRAL APPROACH. 

In this Section we want to study the conductance os- 
cillations by considering the probability amplitude for 
transmission and reflection of electrons and holes between 
the reservoirs via an S-N-S system the diffusive transport 
regime (see Fig. [[]) as a sum of Feynman paths£3 As we 
will show below, one does not actually need to do any 
complicated summations to find this probability ampli- 
tude, because for electron energies below the Thouless 
energy E c (or equivalently for temperatures below the 
Thouless temperature T c = Ec/ks) the hole exactly re- 
traces the electron diffusive path after Andreev reflection. 
It follows that the phase gain of the electron along any 
resonant path between the N-S boundaries (see Fig. |l|) 
is compensated by the hole phase gain along the same 
path. Therefore the phase gain is determined only by 
the phases imposed on the quasiparticles by the super- 
conductors when a trajectory encounters the N-S bound- 
aries. As a result the amplitude does not depend on 
either the form, the length of the diffusive path between 
the superconductors or the configuration of impurities 
(which means there is no need to perform any ensem- 
ble averaging of the conductance). The dependence of 
the resonant probability amplitude on the phase differ- 
ence between the superconductors and on the scattering 
amplitudes at the barriers is easily found by calculat- 
ing the number of reflections at the N-S boundaries and 
the number of backscattering events at the barriers. The 
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conductance is equal to the probability of transmission 
(the modulus squared of the probability amplitude) mul- 
tiplied by the number of different classical resonant paths 
(more strictly, by the number of tubes of width ~ Xp 
around these pathsc3) starting out from a reservoir lead; 
a number that can be straightforwardly estimated. We 
emphasize again that since the conductance associated 
with resonant transmission and reflection does not de- 
pend on the impurity configuration there is no need to 
average it with respect to the impurity positions. 

We start by deriving an equation that connects the 
Feynman path integrals for electrons and holes. To do 
this we shall need the boundary conditions at an N-S 
boundary for the relevant Green's functions. 

The probability amplitude K^'^ (r, r'; t — t') for an 
electron (hole) to propagate from point r at time t to 
point r' at time t' is given by the time dependent Green's 
function satisfying the following equation 

(rihd/dt + H) K {e ' h) (r, r';t-t') = 5{r - r')5{t - t') . 

(59) 

Here the plus (minus) sign is for electrons (holes). The 
initial condition is 



K {e ' h) (r, r'; * - t') = for t - t' < 



(60) 



andi? in Eq. ([39]) is the Hamiltonian describing a metal 
in the diffusive transport regime: 

H = -(n/2m)V 2 + V irnp (r)-e F . (61) 

The potential Vi mp is 

IW(r) = 5>(r - Rj), (62) 

i 

and v(r — Rj) is the potential of an impurity at point 
R 3 . 

In order to derive the boundary conditions we observe 
that the time Fourier transform of K^ e ' h ^ (r, r'; t — t') for 
the case of electrons satisfies the equation 

(H — E — iT))K$ (r, r') = 5(r - r') 



while for the hole case the equation is 

'(h). 



(H + E + iT])K^'{r, r') = S(r - r') 



(64) 



(rj is a small positive constant). At the N-S boundaries 

the Green's functions K E (r, r') are connected with 
each other by the Andreev reflection condition for a fixed 
energy: 

K W E {v^ 2 \v') = e^ + ^^K^\v^ 2 \v') . (65) 

Here r^ 1,2 ) and $12 are the coordinate and the phase 
of the gap function at the first (second) N-S boundary, 



respectivly, and e** E = \A\/(E - i^\A\ 2 - E 2 ), where 
I A I is the magnitude of the gap. Now, an inverse Fourier 
transformation of Eq. ( |65| ) results in the following rela- 
tion, 

dT'i4 e) (r (1 ' 2 V;r') 

-0 

/>oo 

x / dEe** E e i( - T '-^ E / n (66) 



where r = t — t' . We are interested in the case when 
the characteristic time of transmission t<i — t\ is of the 
order of the time, L 2 /D, it takes to diffuse the length L 
of the sample. Since this time can be expressed in terms 
of the Thouless energy as %/E c , the characteristic time 



difference \t — t'\ in the last integral, of Eq. (66) is of the 



order of L 2 /D = h/E c . ThereforeEl in the last integral 
of Eq. d66| ) the main contribution is from energies inside 
an energy interval of order E c <C |A|. In this interval 
E Rj 7r/2 and hence the second integral in Eq. ( pff ) is a 
Dirac (5-function. Therefore the boundary condition for 
the time-dependent electron- and hole Green's functions 
at the N-S boundaries is 

A-W( r (1.2) >r ' ;T ) = ^(0^/2)^)^(1,2)^,.^ (6?) 

According to the Feynman approachB a probability 
amplitude can be written as a path integral. Here we are 
interested in the following probability amplitude, 



K(r 1 ,t 1 ;r 2 ,t 2 ) 



(r 2 ,t 2 ) 



(rut!) 



=l S{r(t)}/fi 



D{r e (t)}D{r h (t)} 
(68) 



This expression sums over all possible paths of an elec- 
tron which start from a point ri in the first reservoir at 
time ti and end at a point r 2 at time t 2 in either the first 
or second reservoir and of all possible hole paths that ap- 
pear due to Andreev reflections at the N-S boundaries. 
For any path the classical action is 



S = / Cdt + *a , (69) 
Jti 

where the phase "J ,4 will be discussed below. The La- 
grangian C for those sections of the path where the par- 
ticle moves as an electron is 



C = mil - V (r e (t)) 



(70) 



For those sections of the path where the particle moves 
as a hole we have 



C = -mr 2 h + V(r h (t)) 



(71) 



Here m is the electron mass, r e and are the elec- 
tron and hole coordinates, respectively, the potential is 
V = Vo + Vi m p] Vq describes the barriers between the 
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sample and the leads to the reservoirs [Vi mp is defined in 
Eq. (|6^)]. While performing the integration in Eq. ( |6^ ) 
one has to use the N-S boundary conditions for electron 
and hole trajectories given by Eq. (|67j). The boundary 
conditions results in an additional term ^,4, which de- 
pends on the macroscopic phases of the superconductors: 



V A = (it/2 + 0!)P« + (tt/2 + fc)P^ 
+ ( 7 r/2- ( ^i 1) + (V2-0 2 )if ) 



(72) 



In this expression -P/H\ and P}^' h \ count how many 



(2) 



(e,h) 



(e,h) 



electron-hole (hole-electron) transformations that has oc- 
curred at the N-S boundaries for a certain trajectory. 

Transport properties of a diffusive system are usually 
calculated in the semiclassical approximation, which im- 
plies (for instance) that the cross section for impurity 
scattering is larger than A^ . We adopt this point of view 
when we now proceed to calculate the functional integral 
in Eq. (68). This means that the method of steepest de- 
scent is useful for performing the integration in Eq. (68), 



and hence classical trajectories thaLjninimize the action 
Eq. ( |69| ) contribute to the integralo 

In order to get the probability amplitude for transmis- 
sion of an electron having fixed energy E from point i"i 
of reservoir 1 to point r 2 of reservoir 1 or 2 at all times 
one has to sum all the relevant amplitudes at all times, 
that is to integrate the amlitude K (multiplied by a fac- 
tor exjp(iEr /Kj) with respect to time r = t% — t%. From 
this we conclude that the probability amplitude A(E) for 
transmission (or reflection) is equal to 



A(E) = J2 R(S)exp 

{S} 



(r 2 ) 



ps{s)ds/h + V A {S) 



(n) 



(73) 



where summation is with respect to classical trajectories 
S that start at point ri and end at point r 2 in reservoirs 
along which an electron with energy E reaches reservoir 
2 (since an electron or a hole) starting from reservoir 1, 
or is reflected back into reservoir 1 (as an electron or a 
hole); ps{s) is the classical momentum as a function of 
the coordinate s along trajectory S, being equal to the 
electron momentum p e and hole momentum p^ at the 
electron- and hole sections of the trajectory S, respec- 
tively. R(S) is a product of the probability amplitudes 
of reflection and transition at the barriers between the 
sample and the reservoirs that occur for the electron and 
hole along path S. ^a(S) is the phase gained along the 
path S by Andreev reflections at the N-S boundaries. 
When counting the number of the trajectories one has to 
take into account the fact that the trajectories has to be 
considered as tubes with a width of the order of the de 
Broglie wavelength A^b = h/pp (see above). 

Now we can calculate the probability amplitudes for an 
electron at the Fermi level with energy E ~ from one 
reservoir to be reflected as a hole back to the same reser- 
voir and to be transmitted as an electron to the other 



reservoir via the diffusive normal metal part of the sam- 
ple. It is crucial for the calculation that at E — under 
Andreev reflection the hole and electron momenta are 
equal but their velocities have equal magnitude but op- 
posite signs. This means that the classical trajectories 
of the electron and hole that end and start at the same 
points at the N-S boundaries exactly repeat each other 
in both ballistic and diffusive samples (as the classical 
trajectory is uniquely determined by the starting point 
and the velocity of the particle). Hence it follows that 
for any classical trajectory with Andreev reflections, at 
E = the total classical action ps(s)ds/h (which is 

the sum of the electron and hole actions) is equal to zero 
as the electron and hole momenta are equal (p e = ph) 
and the integrations are along the same trajectory but 
in opposite directions. Therefore the phase gain along 
such trajectories [see Eq. ( |7^) ] does not depend on ei- 
ther their form, the length of the diffusion path or on the 
configuration of impurities. For resonant transmission 
the summation in Eq. (|7^) with respect to the ampen- 
condlitude scatterings at the junctions is easily carried 
out in the case of low transparency of the barriers at the 
junctions. 

The conductance of a hybrid sample containing both 
normal metal and superconductors, the normal conduc- 
tance is determinecH by the Landauer-Lambgrt formula 
(a), which for a symmetric system reduces toQ 



G 



2e 2 



{T + Ra) 



(74) 



The probabilities Tq and Ra were defined above [cf. 
Eq. (^)J. It is important to note that trajectories which 
connect the two reservoirs, necessarily have a different 
number of electron- and hole sections, while for trajec- 
tories which start and end in the same reservoir these 
numbers are equal. This is a crucial circumstance when 
one sums amplitudes in order to get the total transmis- 
sion amplitude and implies that there is no complete 
compensation of the electron and hole phase gains along 
those trajectories which contribute to the transmission 
amplitude of quasiparticles. As a result destructive inter- 
ference suppresses the transmission amplitude, and the 
main contribution to the conductance is from those tra- 
jectories along which the electron is reflected back into 
the same reservoir as a hole. This is the channel to be 
discussed below. A classical path corresponding to this 
type of reflection is shown in Fig. 0. 

After passing the beam splitter at junction A (that is 
after tunneling through the barrier of this junction) the 
classical diffusive electron trajectory can first encounter 
either the left N-S boundary (clockwise motion) or the 
right N-S boundary (counter-clockwise motion). Adding 
the amplitudes of clockwise and counter-clockwise tra- 
jectories (they form a geometric series) and expanding 
the amplitudes in S(f> = 4> — t <C 1 one finds the total 
probability for an electron being reflected back into the 
same reservoir as a hole to be 
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( 7 2 + <50 2 ) 2 



Here 7 <C 1 is the probability to pass through the barrier 
at the junction. 

From Eq. (75) it follows that the electron-hole 
backscattering amplitude is zero if 4> = tt. This is due 
to the interference between the clockwise and counter- 
clockwise quasiparticle trajectories (in the sense dis- 
cussed above) and can be explained as follows. The am- 
plitude of the electron-hole backscattering can be repre- 
sented as a sum of contributions arising from trajectories 
with different number of Andreev reflections at the super- 
conductors. The ratio between successive terms in this 
geometric series is equal to the amplitude of one Andreev 
reflection at each of the two N-S boundaries. Therefore 
it depends on the phase difference between the supercon- 
ductors and becomes equal to one at resonance, when 
trajectories with very large number of Andreev reflec- 
tions give the same contribution as the ones containing 
only few Andreev events. This is of course the reason why 
a resonance in the electron-hole backscattering channel 
occurs. In addition all terms in the series will be multi- 
plied by a factor e 1 ^ 3 where s labels the N-S boundary 
from which the electron first is Andreev reflected. In 
our notation s — 1 for clockwise trajectories and s = 2 
for counter-clockwise trajectories. An important conse- 
quence of the existence of these multipliers is that on 
resonance, when (4>i — <f>2 = tt), the ratio of this extra ex- 
ponents for clockwise and counterclockwise trajectories is 
equal to —1 and the resonant contributions from clock- 
wise and counterclockwise trajectories to the amplitude 
for electron-hole backscattering cancel each other. A vis- 
ible manifestation of this cancellation effect is a splitting 
of the resonant, conductance peak near <j) — tt so that 
G(cj) = tt) = OB 

If there are Schottky barriers at the N-S boundaries 
additional multipliers appear in the amplitudes for clock- 
wise and counterclockwise trajectories. These are r]l'e l ^ s 
(*) 

S boundary (s = 1,2)). 

In the case of non-equivalent barriers, r A -/- 1 . 
is no compensation of the clockwise and counterclockwise 



(r^ 1 is the amplitude of Andreev reflection at the s:th N- 



(1) ^ A 2) , there 



contributions as is the case when 



— r^?—l. In fact, if 



<C r^' the splitting of the resonant peak disappears. 
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FIG. 12. The resonant conductance peaks are split due to 
interference between clockwise and counter-clockwise quasi- 
particle motion along the trajectories that are associated with 
the formation of Andreev levels. Cf. Eq. (76); here 7 is 0.05 
and 0.2. 

In the semiclassical approximation the total number of 
electrons that contribute to the resonant phase-sensitive 
conductance is equal to the number N± of semiclassical 
tubes of diameter that cover the cross-section of a 
lead between the reservoir and the diffusive sample (as- 
suming the lead has a smaller cross section than the N-S 
boundaries). Hence from the Lambert formula ( ff4| ) and 
Eq. (|75j ) it follows that in the semiclassical approximation 
the phase-sensitive resonant conductance for a diffusive 
sample is equal to 



G = {2e 2 /h)N ± 



7 



(7 2 + H 2 ) 



2\2 



(76) 



Equation ( [76[) implies that the resonant conductance 
peaks are split in such a way that the conductance goe s 
to zero when is an odd multiple of tt (see Fig. fl2| ). 
This splitting appears due to the interference between the 
clockwise and counte-clockwise motion of the particles 
inside the normal sample when electrons are reflected as 
holes back to the same reservoir (see above). 

The above calculations give a qualitative explanation 
to analytical and numerical results for the diffusive case 
presented in Refs. p^|Jl^Jl9| ^2| , ^4| if the results are ob- 
tained for low barrier transparency of the junctions be- 
tween the sample and the reservoirs. 

It should be noted that for the geometry considered in 
some of the papers cited above, where there is only one 
reservoir present, the conductance is determined only by 
the probability for an electron to be scattered back into 
the reservoir as a hole. Therefore the conductance is de- 
termined by the same equation ( |76| ) and hence must also 
be equal to zero at cf> equal to odd multiples of tt for 
the equivalentN-S barrier case. The splitting must dis- 
appear for nonequivalent barriers at the N-S boundaries 
(see, e.g., Fig. 6 and Fig. 7 in Ref. E4; a decrease of the 
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barrier transparency at the junction between the sample 
and the reservoir results in the conductance peaks being 
close to those shown in Fig. ( |T^ ) of this paper). 

We conclude this Section by using the Feynman path 
integral approach to qualitatively consider the temper- 
ature dependence of the oscillating part of the conduc- 
tance for the diffusive case and for temperatures above 
the Thouless temperature E c /ks- If the energy of an 
electron-hole excitation is not equal to zero there is no ex- 
act compensation of the phases gained along the electron- 
and hole portions of the paths connected by Andreev 
reflections. In this case the phase of the transmission 
amplitude A depends on the lengths of the electron-hole 
paths, which in their turn depend on the starting points 
inside the lead between sample and reservoir. The con- 
ductance of the system is a sum of absolute squares of 
amplitudes corresponding to trajectories with different 
starting points in the lead for the classical paths. On the 
other hand classical paths starting from points separated 
by a distance greater than \p, meet different random 
sets of impuritiesEj. As a result their path lengths have 
random values. Hence it follows that one can change 
the summation over starting points to a summation over 
path lengths when calculating the averaged conductance. 
For this purpose we assume a Gaussian length distribu- 
tion of the diffusive paths that start at one N-S bound- 
ary and end up at the other N-S boundary (we choose the 
Gaussian form of the distribution function as an example; 
simple considerations show that choosing a more general 
distribution function only results in additional factors of 
order unity), 

Here Cq is the average length of the paths (Co — 
L 2 vf/D) ■ By averaging the conductance at a fixed en- 
ergy over random path lengths described by the distribu- 
tion function ( J77| ) one easily finds a cut-off factor of order 
exp (—[E/E c f/A\ appearing in the interference terms of 
the conductance. Therefore, destructive interference sets 
in at E ^S> E c (this well known fact justifies the form of 
the distribution function assumed above). The conduc- 
tance oscillations caused by interference only occur for 
energies below or of the order of E c . As a result, at tem- 
peratures T 3> E c /kB the amplitude of the conductance 
oscillations decreases with the temperature as E c /T in 
agreement with Refs. ||]||| 

In order to find the temperature dependence of the gi- 
ant conductance oscillations discussed here when T ^> 
E c /k,B arid for low junction barrier transparencies we 
sum over the paths contributing to the resonance effect 
at a certain energy E, average the conductance over the 
path lengthsusing the distribution function ([77]) and inte- 
grate over energy taking the factor (—dfo/dE) properly 
into account (/q is the Fermi function). As a result we 
find that the oscillating part of the conductance caused 
by the resonant effect is 



SGre. = N ± j^-£g(<l>) (78) 

Here e r is the transparency of the barrier at the junction, 
g(4>) is a 27r-periodic temperature independent function 
with an amplitude of order unity. 

The physical reason for the result (f78|) is that the po- 
sition of the resonant energy peak is tuned by the su- 
perconducting phase difference <j>. With a change of <p 
it can be inside or outside the energy interval of order 
E c associated with the conductance oscillations. As the 
width of the resonant peak is SE ~ e r E c the main con- 
tribution to the conductance oscillations comes from the 
energy interval E ~ t r E Cl and hence the relative num- 
ber of quasuparticles contributing to the oscillations is 
e r E c /ksT. This is why this factor appears in Eq. (f78|). 

VII. CONCLUSIONS 

In this paper we have presented a more thorough dis- 
cussion than in a previous short communicationtS of gi- 
ant conductance oscillations in hybrid mesoscopic sys- 
tems of the Andreev interferometer type, i.e. S-N-S struc- 
tures where the N-part is connected to normal electron 
reservoirs. In Ref. |l8| giant conductance oscillations were 
predicted for a ballistic normal sample when transverse 
mode mixing was absent. The origin of this effect is a de- 
generacy ("bunching" effect) of the Andreev energy lev- 
els at the Fermi energy. This degeneracy of the Andreev 
spectrum arises due to an equality of the longitudinal 
momenta of Fermi energy electrons and -holes undergo- 
ing Andreev reflection. Any process that violates this 
equality lifts the degeneracy and, therefore, decreases the 
amplitude of the conductance oscillations. 

In this paper we considered the effect of giant conduc- 
tance oscillations taking into account transverse mode 
mixing at the junctions between the normal part of the 
sample and the reservoirs. We also considered normal 
reflection in addition to Andreev reflection at the N-S 
boundaries, and scattering of electrons and holes by im- 
purities inside the normal sample. 

Normal reflection of quasiparticles at N-S boundaries 
decreases the propability of Andreev reflection and as a 
consequence also decreases the amplitude of the conduc- 
tance oscillations. We have shown that the probability 
amplitude for the oscillations is giant (that is propor- 
tional to the number of transverse modes N±) until the 
amplitude of the normal reflection is smaller or of the 
same order as the transparency |e r | of the barriers at the 
junctions. 

We have also shown that giant oscillations survive in 
a diffusive sample at temperatures much lower than the 
Thouless temperature. This is because after the electron- 
hole transformation associated with an Andreev reflec- 
tion the electron and the hole move along the same clas- 
sical diffusive trajectory in opposite directions but with 
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equal momenta. As a result the phase gain of the elec- 
tron and the hole along this diffusive path compensate 
each other. The probability amplitude for transmission 
through the sample does not depend on the form or the 
length of the diffusive path, but only on the phase dif- 
ference between the superconductors (i.e., there is no de- 
structive interference). The number of all possible dif- 
ferent semiclassical paths is S/X% — N±, where S is the 
cross-section area, as each path has a width of the order 
of the de Broglie wave length Xb- Therefore the ampli- 
tude of the conductance oscillations in the diffusive case 
remains giant and proportional to the number of trans- 
verse modes N± as for ballistic samples. The above qual- 
itative picture agrees with analytical calculations for the 
diffusive case by Beenakkei; et al&3, Zaitsevliil, Alkopp et 
alxl, Volkov and Zaitsevu, and Claughton et aZ.E3. 

The presence of potential barriers at the junctions be- 
tween the sample and the normal electron reservoirs is 
most crucial for the giant oscillations to exist. In the 
weak coupling case (e r <C 1) we have shown analyt- 
ically and numerically that the amplitude of the con- 
ductance oscillations is independent of the barrier trans- 
parency and proportional to the number of transverse 
modes Nj_ . When the transparency of the barriers is in- 
creased, our numerical calculations show that the ampli- 
tude decreases. In the absence of barriers at the junctions 
the amplitude becomes zero. The latter result agrees 
with the sum rule in Ref. ^l] according to which the con- 
ductance is equal to the number of transverse modes — 
in the absence of barriers — times the conductance quan- 
tum and does not depend on the superconducting phase 
difference. This can be qualitatively understood since an 
electron (hole) coming from the reservoir after being first 
Andreev reflected at one N-S boundary as a hole (elec- 
tron) then at zero temperature returns to the reservoir 
by retracing the path of the incoming particle without 
reaching the second N-S boundary. n 

Recently Nazarov and Stoofl, and Volkov et alB pro- 
posed a new mechanism for conductance oscillations in 
diffusive samples that is effective if the temperature is 
close to the Thouless temperature (thermal effect). They 
used the dependence of the diffusion coefficient on the 
quasiparticle energy and found conductance oscillations 
with the superconducting phase difference in the absence 
of barriers. The amplitude of the oscillations was found 
to reach its highest value at the Thouless energy. We 
propose that this effect can be qualitatively understood 
if one takes into account the fact that under Andreev 
transformation at an N-S boundary there is a finite an- 
gle between the trajectories of the incident- and reflected 
particle, which is proportional to their excitation energy. 
Simple estimations show that near the Thouless tempera- 
ture the classical trajectories of the electron and the hole 
can be separated by a distance of the order of the de 
Broglie wavelength (that is the width of the semiclassical 
trajectories) before the particle leaves the normal diffu- 
sive part of the sample for a reservoir. When this happens 
the trajectories meet different sets of impurities and can 



diffuse along very different paths inside the sample. This 
permits the quasi-particles to encounter also the other 
N-S boundary and undergo Andreev reflection there be- 
fore leaving the sample. Therefore the conductance starts 
to depend on the superconducting phase difference and 
conductance oscillations arise. When the temperature is 
higher than the Thouless temperature, the phase gains 
of the quasi-particles along the semiclassical paths are 
not compensated and are much larger than unity; the 
destructive interference kills the thermal effect, in agree- 
ment with the results of the papers cited above. 

When the transparency of the barriers at the junc- 
tions has intermediate values both the thermal effect and 
the resonant oscillation effect considered in this paper 
are in effect simultaneously provided the temperature is 
near the Thouless temperature. The effects can be dis- 
tinguished by decreasing the temperature, which results 
in a decrease of the amplitude of the conductance os- 
cillations in the case of the thermal effect (vanishing at 
zero temperature) while the resonant amplitude of the 
conductance increases and has its largest value at zero 
temperature. 

This work was supported by the Swedish Royal 
Academy of Sciences (KVA) and by the Swedish Natural 
Science Research Council (NFR). 

APPENDIX 

A. Elements of the scattering matrix describing 
coupling 

We are to describe leads of finite width and the scat- 
tering matrix should mix different modes. The S matrix 
will then be of size 3A^j_ x 3N± and the unitary condition 
for the submatrices of size N± x N± gives 



sns T n = 1 - 2s 12 s T 12 (79) 

S22S22 + S23S23 = 1 - si2s\ 2 (80) 

S22S23 + S23S22 = -h2s\ 2 (81) 

*ii«Ia + »u (4 2 + 4») = ( 82 ) 

S12S11 + (S22 + S23) s\ 2 = (83) 



Hence we have 5 matrix equations for 8 matrices (as each 
matrix has an independent hermitian and antihermitian 
part) there are 3 undetermined matrices. We choose 
them to be §12 and the antihermitian part of §22- 

As it was said above the antihermitian part is not de- 
termined by the unitary conditions for the matrix S. 
In general case it is of the same order of magnitude 
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as the hermitian part as they are connected with the 
Kramers-Kronig relation. Therefore the matrix elements 
of S22 oc e r . An analogous analysis of the rest of the 
equations shows the matrix elements of §23 — 1, Pu and 
V to be of the same order of magnitude. 

J A) 

If S12 and the anti-hcrmitian part s 22 commutes and 
other matrices are expressed in terms of §12 they may 
be simultaneously diagonalized. The N± eigenvalues of 
matrices Sij are denoted Ay were indices numbering the 
eigenvalues are suppressed. Directly from the unitary 
conditions we get 



|Aii| = A /l-2|Ai 2 | s 



IA23I = v /1 -| A i2l 2 -|A 2 2| 2 
A*iAi2 + A^ 2 (A22 + A23) = 



|A 12 | 2 + 2|A 22 ||A23| cos I 



»22 



•'23 



) = o 



(84) 
(85) 
(86) 
(87) 



which gives the requirement 022 — ^23 € [^/^i Stt/2]. Now 
use a hermitian §12 i.e. real eigenvalues, 0x2 = Tin. Put 
023 = in order to agree with the weak coupling limit 
where no phase gain is expected in passing the reservoir. 

(A) 

We get with A12 and A 22 as eigenvalues of known her- 
mitian matrices 



|Au| = Vl-2|Ai 2 | 2 

.(A) _ ,(A) 
A ll — A 22 

(H) _ 



X 



11 



Ml — A 11 + tA xl 



v(H) 
^22 



A22 — A 22 + «A 22 




A23 = a/i- IA12P- |A 22 | 2 



(89) 
(90) 
(91) 

(92) 

(93) 
(94) 



A symmetric random matrix with normally distributed 
elements will have real eigenvalues distributed according 
to the semicircle law. We have used values of mean and 



(A) 



variance 1 creating random matrices S12 and s 22 

Sjj = (4>n\Sij\4> m } 

where \4> m ) is a complete set of vectors 

N ± 
a=l 



(95) 



(96) 



where \tp a ) are eigenvectors to Sy. If the corresponding 
eigenvalues are called A a 



N± N± 
Jnm) _ \ ^ \ ^ (to) An) 



N± 



a=l /3=1 



With elements in the matrix §ij of order unity 



(97) 



)|2 _ 



(98) 



„( m ) 



we get 7q" v cx l/\/N±. Our random matrix S12 will 
be multiplied by \J e,./iVj_ before eigenvalues are calcu- 
lated. Then by setting I = e r /e c = 1 and approaching 
the strong coupling limit a maximum value e c is found 
fulfilling the unitary conditions. The strength of the cou- 
pling of the reservoirs is now parametrized by e S [0,1]. 
The matrix S22 is in the weak couoling limit for one chan- 
nel seen to be proportional to e r c3 and therefore the ran- 
dom matrix giving A 2 2 i s multiplied by e r /y/N±. Then 
by using the eigenvectors of the matrix s\2 we transform 
all the matrices Sy back to the initial representation in 
which they are not diagonal and their matrix elements arc 
the probability amplitudes of scattering to the respective 
transverse modes. 

To realize another type of scattering matrix to describe 
the coupling to the reservoirs we do as follows. The 
essential features of the junction are coupling to elec- 
tron reservoirs and mixing between modes, both features 
may be parametrized by the strength of the coupling e. 
The coupling to the reservoirs is described by Biittiker 
matrices.Cj If this is done mode by mode there will be 
no mixing. An additional unitary matrix is used to mix 
modes. This matrix has all diagonal elements equal to 
each other and all off-diagonal elements equal to each 
other, describing scattering into the same mode and mix- 
ing between modes respectively. By keeping the elements 
equal an isotropic situation is simulated where scattering 
into any mode is possible. The elements u of the N± x N± 
unitary matrix must fulfill 

\u H \ 2 + (N x - l)\ Uij \ 2 = 1 (99) 

uJitiy + + (M - 2)|«ij | 2 = (100) 

This gives \uij\ 2 < 1/(N± — 1), the phases of Uij and ua 
are connected by 



(JVx-2) K | 



(101) 



we note that for N± > 4, \uij\ and \ua\ are not allowed 
to have the same value. We wish to consider large N_\_ 
and write the elements 



\uu\ = 1 - ce 



N ± - 1 



(102) 
(103) 



in order to agree with the limit e = when mixing is 
expected to be absent since the waves in the decoupled 
sample do not feel the reservoirs. The condition Eq. (Uni- 
tary U) gives 



4(ATj_ - 1) 



(104) 
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The parameters are the coupling e, the number of modes 
N± and the angle <pa. We write the angle <pa £ e[0, 2ir] to 
agree with the weak coupling limit where no phase gain 
is expected. By using different angles we get an ensemble 
of matrices describing different samples. 

The eigenvalues of the unitary U all have amplitudes 
of length unity. This means that all modes will be open 
for transmission .EH Opening ofj-transmission channels has 
been observed in experiments Ej ._. 

To describe coupling the results by Buttikenfj are used 
in diagonal matrices. All these matrices Sy are multiplied 
by U. 



C. Active channels 

As for N±_ 1 the set of equations Eq. (|2l]) can not 
be analytically solved and the amplitudes j n can not be 
explicitly found we estimate the number Nr of transverse 
modes inside the resonant region (—e r hvF / L, e r hvF / L) 
and use Eq. ( ^9|) to get the conductance to within a factor 
of the order of unity. We determine Nr in the following 
way. 

oc 

N «= J dE (EL/ni)i + ef W < 112 > 



B. A barrier at the NS-interface 

A Schottky barrierll at the interface of a semiconduc- 
tor and a metal can be characterized with a scattering 
matrix 



a(o) _ 



ro to 
~to r 



M 2 + |t | 2 = i 



(105) 



that connects the constant factors a\ and b\ of the plane 
waves coming in and going out of the barrier inside the 
semiconductor, respectively, with those a-i and 62 inside 
the metal: 



b i = £4%, 
k=l 



i = 1,2 



(106) 



(Hence \to\ 2 is the transparency of the barrier) When 
the metal is in the superconducting state the matrix of 
reflection of the semiconductor charge carriers at the N- 
S boundary (the semiconductor is on the right and the 
superconductor is on the left of the N-S boundary) is 



rN r A 

~ r A r *N 



where 



r 4 = it f . 



|io| 4 + 4|r | 2 sin 2 ^ E 



e iif>E - | ro |2 e -«A f 



2 sin %l>i 



y / M 4 + 4|r | 2 sin 2 i> E 

e i<t> 

^/M 4 +4|r | 2 sin 2 xj, E 



E - iy/\A\ 2 - E 2 



(107) 



(108) 



(109) 



(110) 



(111) 



Here |A| and 4> are the modulus and the phase of the su- 
perconducting gap, respectively, E is the electron energy 
measured from the Fermi level. From Eq. (107-111) it is 
straightforward to see that \ta\ 



\tnV = 1. 



where v(E) is the state density function 

N± N ± 



I n=l 



, dQ n 
1 dE 



\S(Q n ); (113) 



Here the spectrum function Q n is determined by Eq. (|44|). 

To find the state density function v(E) we use the 
method developed in Ref. |52|. As = ±(hv)~ 1 L with 



v = Ti^kp — k±(n) 2 /m the factor 8Q n /dE in ( |113| ) is a 
trigonometrical function of ±<p± as well as Q n is, and it 
is productive to expand v into Fourier series in (p± and 
write it as follows. 



N± 



n—0 S— — 00 k- 



A 



s.k 



(2tt) 



27T p2-K 

d(p + I dtp- 

u JO 



,9Qn 

1 dE 



\S(Q n )e 



(114) 



-i(stp~+k<p+) 



(115) 



In this paper we assume the length L and the width d of 
the sample to satisfy the inequality 



L 



> 



(116) 



Usin g th is inequality, Eq. (122) and the estimation of 
Eq. (127) in Appendix D one sees that th e ma in contri- 
bution to the state density function Eq. (113) is of the 
terms 



A s .o = 



L 



(2w) 2 h 
x<5(cos (p 



- dip+ (^_|sm^_|e-^- (117) 
" Jo Jo 



\tni \\rN2 I cos lfi + + \tai \\rA2\ COS ( 



Performing inte gration with respect to <3 _ in ( j l 17[) and 
over E in (112) with application of (114), jllq ) one ob- 
tains the conductance 



E 



,-2|a|e 



cos sip\((p+)d(p+ (118) 
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ipi{(p+) = arccos(jrjvi||rjv2l C0S ( i P+) ~ k^ill^l cos^J 
if < tp_ < 7r, and 

<pl(<p+) = 27r-arccos^|r A r 1 ||r A r 2 | cos(^ + )-|r yU ||r j42 | cos</>) 



if 7r < < 2it. Performing summation in (118) one gets 

e 2 

G = Nj_ — e r x 
irn 



(119) 



2ts 



e r dif + 



o 1 - |nvi||rjv 2 | cos</3+ + |r J 4 1 ||T-A 2 | cos0 + e 2 /2 



as l r iVx 2I < 1 integration in (119) gives 



therefore, S is less than N_\_. As Nj_ ~ kpd we have 
the following estimation of the sum of fast oscillating 
terms. 



S < y/Lk^/(k F d) 



(127) 



Therefore the sum of fast oscillating terms can be ne- 
glected if 



> 



(128) 



G = N 1 



7r?i^(l + |r Al ||r A2 |cos</) + e2/2)2-| rAri |2| rjV2 |2 ) 

(120) 



If the boundaries are symmetric, that is rjvi = r/V2 7 the 
conductance is 



G = N 1 



^ V / (1 + M 2 COS0 + £ 2/ 2 )2 _| rjv |4) 



(121) 



D. Fast oscillating terms 

In this Appendix we evaluate a sum of fast oscillating 
terms 



s= -L^r e ^i~(^) 2 ] 1/2 



(122) 



that appears in the density state function Eq. (113) and 
the transition probability Eq. (|27|), Eq. (|2|). Using the 
Poisson formula one can write 



s = urn — y [ dxe ^v^y^+^ e -n\ k \s 



k= — oo 



(123) 



Here 



X = Lk Fl a = k F d/ir (124) 
As A ^> 1 one can use the saddle point method to get 

(125) 

A 2 



5 = -L^v^-" /4 x 

A ±% /A 



00 / X 2 \ 3/4 / 



k— — OQ 



From Eq. (125) it follows that terms with k which are 
less or of the order of A/a contribute to the sum and, 
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